Spatially resolved spectroscopy of monolayer graphene on Si02 
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We have carried out scanning tunnehng spectroscopy measurements on exfohated monolayer 
graphene on Si02 to probe the correlation between its electronic and structural properties. Maps 
of the local density of states are characterized by electron and hole puddles that arise due to long 
range intravalley scattering from intrinsic ripples in graphene and random charged impurities. At 
low energy, we observe short range intervalley scattering which we attribute to lattice defects. Our 
results demonstrate that the electronic properties of graphene are influenced by intrinsic ripples, 
defects and the underlying Si02 substrate. 
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Graphene is the two dimensional form of carbon char- 
acterized by a honeycomb lattice with two inequiva- 
lent lattice sites. It is unusual in many aspects, one 
of them being the linear energy dispersion relation that 
causes the electrons to obey the relativistic Dirac equa- 
tion instead of the Schrodinger equation^. The theo- 
retical prediction of the instability of 2D crystals along 
with unsuccessful efforts at their synthesis on insulat- 
ing substrates kept graphene far from the realm of ex- 
periments. A breakthrough in the isolation of graphene 
from 3D graphite opened up a new frontier for exper- 
imental investigation^'^. Electrical transport measure- 
ments have brought to the fore some of the exotic prop- 
erties of graphene like a novel quantum Hall effect, high 
carrier mobility and minimum conductivity^i^i^. Trans- 
mission electron microscopy characterization (TEM) of 
graphene revealed intrinsic ripples^ which were confirmed 
by Monte Carlo simulations^. These structural findings 
hinted that electron scattering and carrier localization 
could be highly influenced not only by external effects 
such as impurities but also by intrinsic ripples. Further- 
more, electrical transport measurements showed the sup- 
pression of electron localization^ and that the mobility 
was insensitive to doping^^. Both these results, unex- 
pected for a two dimensional system like graphene, were 
attributed to corrugations in graphene. Hence, to in- 
vestigate phenomena involving electronic and structural 
aspects of graphene we have performed spatially resolved 
spectroscopy measurements. 

A single atom thick graphene sheet can be studied by 
suspending it on a micron-sized metal grid or by support- 
ing it on a substrate like SiOsr^iiiii^. The sheet is not 
perfectly flat but contoured with intrinsic ripples with 
deformations up to 1 nm normal to the plane clearly evi- 
dent from TEM of suspended graphene and Monte Carlo 
simulationsii^. These ripples are due to the stability re- 
quirement of the two dimensional lattice and the ability 
of carbon to bond with asymmetric bond lengths^'^. Sup- 
ported sheets of graphene have an additional constraint 
on their morphology due to the Si02 substrate. They 
have been found to partially conform to the substrate 
while still having intrinsic ripples and therefore areas not 



in contact with ^iO^^-^^. These previous measure- 
ments showed the topographic structure of monolayer 
graphene. Spectroscopy measurements of graphene on 
Si02 revealed the contribution of phonons to the tun- 
neling processi^. However, none of these measurements 
probed the effects of graphene morphology on its elec- 
tronic properties. Theoretical calculations show that 
variations in the local curvature of graphene, ripples, lead 
to a TT — <j orbital mixing which changes the local electro- 
chemical potential breaking particle hole symmetry and 
causing charge inhomogeneities^^. The shift in the local 
chemical potential $(r) is proportional to the square of 
the curvature of graphene 
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where h{r) is the local height referenced to a flat con- 
figuration, a is the nearest neighbor distance and a ^ 
9.2 eV is a constanti^. This implies that curvature of 
the graphene sheet leads to density variations due to the 
shifting chemical potential. 

In addition to the intrinsic rippling in graphene, its 
two dimensional nature makes its entire surface area sus- 
ceptible to adsorption and defects. Given the intricacies 
involved in graphene device preparation the presence of 
random impurities, adsorbates or defects from various 
sources is difficult to prevent. Depending on the type of 
imperfection, they can either act as long or short range 
scattering sites. Coulomb scattering by random charged 
impurities at the graphene-Si02 boundary is one exam- 
ple of long range scattering. This long range, intravalley, 
scattering tends to create inhomogeneities in the carrier 
density forming electron and hole puddles^^'^^. When 
the scattering potential has components shorter than the 
lattice constant, typically due to lattice defects, interval- 
ley scattering takes place which mixes the two sublattices 
of graphene^^. Here we use a low temperature scanning 
tunneling microscope (STM) to probe the morphology, 
local electronic properties and scattering phenomena in 
monolayer graphene. 

Graphene was prepared by mechanical exfoliation of 
graphite on Si02^s^. Monolayer areas were identified us- 
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FIG. 1: (a) Schematic of a mechanically exfoliated monolayer 
graphene flake (gray) with gold electrodes on a Si02 sub- 
strate, (b) STM topographic image of monolayer graphene 
showing the hexagonal lattice as well as the underlying sur- 
face corrugations. The scale bar is 2 nm. The imaging pa- 
rameters are sample voltage Vs = 0.25 V, tunneling current 
I, = 100 pA. 



ing an optical microscope and then Ti/Au electrodes were 
deposited using standard electron beam lithography. The 
lithography process leaves some PMMA resist on the sur- 
face of graphene. To eliminate the resist the device was 
annealed in argon and hydrogen at 400 °C for 1 hour- 
followed by annealing in air at 300 °C for 30 minutes. 
The device was then immediately transferred to the STM 
(Omicron low temperature STM operating at T = 4.5 K 
in ultrahigh vacuum (p < 10~^^ mbar)). Electrochem- 
ically etched tungsten tips were used for imaging and 
spectroscopy. All of the tips used were first checked on 
an Au surface to ensure that their density of states was 
constant. 

A schematic of the measurement set-up showing the 
graphene flake on Si02 with a gold electrode for electri- 
cal contact is shown in Fig. [TJa). A typical STM image 
of the monolayer graphene showing the complete hexag- 
onal lattice along with surface corrugations due to the 
underlying Si02 substrate is shown in Fig. [TJ^b). The ob- 
servation of both sublattices, giving the hexagonal struc- 
ture, is characteristic of monolayer graphene as opposed 
to bulk graphite where only a single sublattice is usu- 
ally observed giving a triangular pattern. The hexagonal 
lattice in our images extends over areas as large as 40 
nm demonstrating the cleanliness of the surface. The 
uneven surface underneath the hexagons is a distinctive 
feature of all our images. This additional corrugation has 
a height variation of ~ 5 A over an area of 30 x 30 nm^. 
The lateral extent of these corrugations is in the range of 
a few nanometers mimicking the Si02 corrugation^^'^^. 

We have recorded the local density of states of 
graphene using dl/dV point spectroscopy measurements. 
In this case, the tip is fixed at a specific location on the 
sample, the feedback loop is turned off and the sample 
voltage is ramped within a specific energy window. An 
ac modulation voltage of 4 mV rms, 574 Hz is applied to 
the sample and the resulting dl/dV spectrum is recorded 
using lockin detection. One such spectrum is shown in 
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FIG. 2: dl/dV point spectroscopy shows a linear relationship 
between the tunneling conductance over a sample voltage Vs 
of dz 0.15 V. The inset shows the region near the Fermi level, 
no gap is seen. 



Fig. [21 The dl/dV spectrum is approximately linear in 
energy and featureless. Neither a gap nor a dip is seen 
around the Fermi level. 

To understand the effect of the corrugations seen in 
Fig. [Hb) on the density of states of graphene, we have 
performed spatially resolved scanning tunneling spec- 
troscopy (STS) measurements. By measuring the dif- 
ferential conductance, dl/dV as a function of energy and 
position, we obtain maps of the local density of states 
(LDOS). Figure [3] shows atomically resolved topography 
and spectroscopy over a 40 x 40 nm^ area. The topogra- 
phy of the area. Fig. [3](a), demonstrates the atomic reso- 
lution as well as the corrugations. Fig. [3Kb)- (e) are maps 
of the LDOS over the same area for four different sample 
voltages ranging from -0.2 V to 0.4 V. The maps show ar- 
eas of high and low differential conductance which change 
as a function of energy. These changes demonstrate that 
the Dirac point varies as a function of position as dis- 
cussed below. At low energy, the shifting of the Dirac 
point leads to electron and hole puddles. 

When performing the spectroscopy measurements, the 
feedback circuit is turned off and the tip height is held 
constant. The height is determined by the setpoint cur- 
rent, /g, (100 pA in this case) and the setpoint voltage, 
Vs^ (0.4 V in this case). If we assume that the Dirac 
point is at the Fermi energy, then the dl/dV curve can 
be written as ^ = where V is the sample voltage 
(green curve in Fig. [4j). The value of a is determined 
by the setpoint parameters. Is and Vg. This is because 
the feedback circuit adjusts the tip height such that the 
integrated current from to Vs is /g. The value of a is 
determined by setting the area under the dl/dV curve 

equal to /g , /g = . This gives a value of a = y^- and 

therefore ^ = ^tt |V| (green curve in Fig. [D^a)). 




FIG. 3: (a) STM topography showing the lattice and sur- 
face corrugations. The scale bar is 8 nm. The topography is 
recorded with Vs = 0.4 V and L - 100 pA. (b) to (e) dl/dV 
maps at sample voltages -0.2 V, 0.0 V, 0.2 V and 0.4 V re- 
spectively showing electron and hole puddles. For all of the 
images, the current was stabilized at the parameters of (a) and 
the feedback was then switched off. The maps were recorded 
using lock-in detection, (f) Spatially averaged dl/dV curves 
at 5 different regions of the graphene indicated in (a). The 
curves show a crossover near 0.2 V as a result of the shifting 
Dirac point. 



However, in a region where the Dirac point shifts by 
an amount A < the sample essentially becomes more 
conductive. Therefore, the tip must move farther from 
the surface to maintain the setpoint current causing the 
dl/dV curve to have a lower slope. Now, ^ = P\V — A\ 
as shown by the blue curve in Fig. [41(a). The value of 
P is once again determined by the requirement that the 
area under the ^ curve from to is /s- This gives 
f3 = y2^2Ay • point where two dl/dV curves with 
slopes a and /3 cross can be obtained by setting the dl/dV 
curves to be equal and solving for V. The result is that 
the two curves intersect at Vs/2. 
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FIG. 4: (a) Illustration of effecting of shifting Dirac point on 
dl/dV curves. The shift in the Dirac point causes the slopes 
of the dl/dV curves to change leading to a crossover at Vs/2. 
Furthermore, the value of dl/dV at Vs gives information about 
the amount of the shift in the Dirac point, A. (b) Standard 
deviation of dl/dV for 3 different locations and values of Vs. 
The setpoint voltage Vs is 0.4 V (blue curve), 0.25 V (green 
curve) and 0.15 V (yellow curve). The minimum occurs at 
Vs/2 for each of the curves. 

When the Dirac point shifts by an amount A > 0, 
^ = j\V — A\ and a similar effect occurs except the tip 
moves closer to the surface and the slope of the dl/dV 
curve increases (red curve in Fig. [H^a)). The value of 7 
is given by 

^ 2Is 
^ - 2AVs + 2A2 

Again the dl/dV curves cross at ^5/2 if A << and at 
negative sample voltages, the curves have different slopes 
based on the energy of the Dirac point but they do not 
cross. 

This analysis makes two predictions for the LDOS 
maps. (1) The variation in the maps should be small- 
est at Vs/2 because regardless of the energy of the Dirac 
point, all the dl/dV curves have the same value at Vs/2. 
(2) The value of dl/dV at Vs is related to the energy 
of the Dirac point by the relationship A = Fg( "^~_^^^ ), 
where r] is the normalized differential conductance given 

u dl/dV 

We have studied the variation in the LDOS maps for a 
set of three different locations and tunneling parameters. 
The results are shown in Fig. HKb). The blue curve is for 
the images shown in Fig. [3] where Vs was 0.4 V. From the 
curve it is clear that the minimum variation occurs at 0.2 
V as can also be seen in Fig. [31(d). The green curve is 
taken in a different region of the sample with Vs = 0.25 
V. The minimum has shifted to lower energy compared to 
the blue curve and it is now at 0.125 V. Lastly, the yellow 
curve has Vs = 0.15 V and the minimum is at the lowest 
energy of the three curves. This shows that the location 
of the minimum in the LDOS maps is dependent on the 
parameters used for the measurement but the value of 
dl/dV at ^5 still gives information on the Dirac point. 

Fig. [31(f) shows the dl/dV curves taken from the LDOS 
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FIG. 5: (a) STM topographic image of 30x30 nm^ area taken 
with Vs = 0.15 V and L = 50 pA. (b) Shift in local electro- 
chemical potential calculated from the curvature of the topo- 
graphic image shown in (a), (c) LDOS map recorded at Vs 
= 0.15 V with the tip height determined using the imaging 
parameters from the topography measurement in (a). Scale 
bar in all images is 3 nm. 



maps for five different areas of the image as shown in Fig. 
[3fa). These curves exhibit a range of dl/dV values at Vs 
which is evidence for the formation of electron and hole 
puddles. In Fig. El^e), the red regions (large dl/dV) 
are locations where the Dirac point has shifted towards 
positive sample voltage while in the blue regions (smaller 
dl/dV) the Dirac point has shifted towards negative sam- 
ple voltage. From the changes in the dl/dV curves we 
estimate that the shift in the Dirac point, AE^^ is about 
77 mV in our images. 

As discussed earlier in this manuscript intrinsic ripples 
in graphene give rise to a spatially varying electrochemi- 
cal potential (equation (1)). This equation predicts that 
highly curved regions of the sample will be electron doped 
while flat regions will be hole doped. Using a STM to- 
pographic image, Fig. [5](a), we have calculated the shift 
in the electrochemical potential caused by the graphene 
curvature using equation (1). The results in Fig. ^h) 
show areas of high curvature, large negative change in 
chemical potential, in blue while relatively flat regions 
are brown. A comparison of the shift in Dirac point from 
the dl/dV map at 0.15 V, Fig. [SJc) and the electrochem- 
ical potential landscape Fig. ^h) shows there are limited 
regions which are in agreement between them. For exam- 
ple the curved region in the bottom right side is electron 
doped (blue region in Fig. [5fc) and Fig. El^b)). How- 
ever, this is not always the case; there are curved regions 
such as the blue patch in the center of Fig. [51(b) that 
is hole doped as seen in the center of Fig. [51(c). Thus 
from our measurements we conclude that the curvature 
in the graphene flake contributes to a variation in the 
electrochemical potential but it is not the main factor re- 
sponsible for the features in the dl/dV map. Instead, the 
potential variation is due to a combination of the ripples 
and long range scatterers such as random charged im- 
purities present on the graphene sheet These electron 
and hole puddles, a signature of disorder in graphene, are 
also responsible for the finite minimum conductivity at 
the Dirac pointii^. 

A second source of potential variation in graphene 
is long range scattering from random charged impuri- 
ties present at the graphene substrate interface. From 



our measured shift in the Dirac point, AE^^ we cal- 
culate the density variation using An = '^^f^a to be 

4.3 X lO^^cm"^. An independent measurement of the im- 
purity density can be done using the voltage dependence 
of the conductivity, a = 20 e eVg / {hn it) where e is the 
dielectric constant of Si02, t is the oxide thickness and 
Vg is the gate voltagei^. We measured another flake of 
graphene processed in the same manner as the one used 
for the STS measurements and found = 4 x lO^^cm"^. 
This impurity density can be converted to local electron 
density fluctuations using Sn'^ = ni/{S7T(P) where d is 
the distance of the impurities from the graphene^^. This 
agreement between our spectroscopy measurements and 
transport measurements as well as transport measure- 
ments by other groupsi^ indicates that the role of charged 
impurities is critical for understanding current graphene 
devices. 

While long range scatterers give rise to intravalley 
scattering, within one sublattice, creating the electron 
and hole puddles seen in the LDOS images, short range 
scatterers such as lattice defects can also be present in 
graphene. These short range scatterers induce intervalley 
scattering from one Dirac cone to the other. The mobil- 
ity of graphene was found to be insensitive to doping 
with gas molecules, long range scatterers^^. This implies 
that there is another source which may already limit the 
mobility of the graphene such as short range scatterers 
which tend to give a constant resistivity and hence low 
mobilit y ^^1^"*^ . We observe lattice defects in the mono- 
layer graphene images and analyze the resulting Fourier 
transforms of the LDOS maps and topographs for an in- 
sight into intervalley scattering. In the case of epitaxial 
graphene and graphite, lattice defects have been shown 
to give rise to scattering and interferenc o^^i^^ . Figure [H 
shows the results for graphene on Si02. Atomically re- 
solved topography. Fig. [Gja), shows the hexagonal lattice 
of graphene and the corresponding Fourier transform. 
Fig. [IKc), contains the reciprocal lattice points (red cir- 
cles) and components due to the C-C bonds (brown cir- 
cles) which are longer and rotated by 30 degrees. Figure 
[61(b) is an image of the LDOS at the Fermi level showing 
the defect induced interference pattern superimposed on 
the electron and hole puddles due to the long range scat- 
tering. The Fourier transform of this image. Fig. [61(d), 
shows two hexagonal patterns. The outer hexagon, red 
circles, is due to the reciprocal lattice and is located at 
the same points as in Fig. [BJc). There is also an addi- 
tional inner hexagon, blue circles, that is rotated by 30 
degrees due to scattering. This is the >/3 x \/3 R30° inter- 
ference pattern. The six peaks at the K± points in Fig. 
[61(d), are a result of intervalley, short range, scattering 
events. 

To characterize the strength of intervalley scattering as 
a function of energy, we have fit the six peaks at the K± 
points of the \/3 x \/3 R30° interference pattern with a 
2D Lorentzian to determine the area under the peak. We 
do not observe any structures around these K± points un- 
like those observed in case of graphene on SiC^^. Hence 
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FIG. 6: (a) STM image of graphene, the scale bar is 6 nm. 
Inset is a 4 nm region within the image showing the hexago- 
nal lattice, (b) LDOS map at the Fermi energy (0 V) taken 
simultaneously. Inset corresponds to the same 4x4 nm^ area 
as in the topography image showing a different pattern due to 
scattering from defects, (c) 2D Fourier transform of (a) show- 
ing the reciprocal lattice (red circles) and components due to 
the C-C bonds (brown circles). The scale bar is 4 nm~^ . (d) 
Fourier transform of (b). The inner hexagon (blue circles) 
represents the interference pattern due to intervalley scatter- 
ing. The outer hexagon (red circles) is due to the reciprocal 
lattice. 



we have chosen a Lorentzian for the fit. This area is 
then divided by the average value of the Fourier trans- 
form to get the relative amount of intervalley scattering 
at a given energy. As an example, Figure [TJa) shows 
one such Fourier transform with six peaks that are fit 
with a 2D Lorentzian. Figure [TJb) shows the linescan 
across one of the 6 peaks in the Fourier transform and 
the Lorentzian fit corresponding to that peak. Figure 
[Tl^c) plots the relative strength of these peaks as a func- 
tion of sample voltage for three different regions of the 



sample. The broad peak in scattering at the Fermi energy 
is independent of the imaging parameters. It is consis- 
tently observed in areas of the flake with scatterers or 
defects. As we move away from the Fermi energy the 
strength of scattering decreases sharply. This enhanced 
intervalley scattering at low energy is evidence of weak 
localization of carrier o^^i^^ . 

In conclusion, we have presented an extensive to- 
pographic and spectroscopic investigation of monolayer 
graphene on Si02 using a STM at 4.5 K. We were able 
to atomically resolve large areas of exfoliated monolayer 
graphene, record energy resolved local density of states 
maps and interpret the maps to be a signature of the 
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FIG. 7: (a) 2D Fourier transform of a LDOS map recorded 
at the Fermi energy. The traces represent the contours cor- 
responding to a 2D Lorentzian fit. The scale bar is 3 nm~^. 

(b) Linescan across one of the peaks in the Fourier transform 
(blue trace) and the corresponding Lorentzian fit (red trace) . 

(c) Amount of intervalley scattering as a function of sample 
voltage. Each curve is from a different region of the graphene 
flake. The blue curve corresponds to the area shown in Fig. [3] 
and is taken with 1^ = 0.4 V. The green curve and the yellow 
curve are from different (30 x 30 nm^) areas taken with Vs — 
0.25 V and 0.15 V respectively. 



shifting Dirac point. Also we identified and analyzed in- 
tervalley scattering mechanisms on exfoliated monolayer 
graphene. Thus, the electronic properties of monolayer 
graphene are closely related to intrinsic ripples, Si02 sub- 
strate morphology and random impurities. Investigation 
of monolayer graphene flakes can be extended to different 
geometries and substrates to advance our understanding 
of the subtleties of the electronic properties of graphene. 
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